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Abstract 

We study the equations of a two dimensional incompressible Newtonian fluid cou- 
pled with a dispersive parabolic-elliptic system on bounded domains. Global in time 
weak solutions are shown to exist and converge with a rate to the stationary solution 
for LP' initial data. This paper extends and improves on a body of work surrounding 
the Debye-Hiickel system to the hydrodyanamical case. 

1 Introduction 



The equations of electro-hydrodynamics are 

du 
'di 



du 

— + u-Vu + Vp = Au + A(j)V(j), (1) 



V • n = 0, (2) 

dv 

— + n • = V • {Vv - vV(j)) , (3) 
dw 

— + u- Vw = V ■ {Vw + wV(j)) , (4) 

A(j) = V — w. (5) 

On bounded domains, solutions are determined by the conditions 

dv d(j) 

^ ~ '"IT = 0' (6 
du ov 

dw dd) 

— + u;/ = 0, 7 

OV Of 

where v is the outward pointing normal to dU and 

n = (8) 



on dU. In this paper solutions of the Poisson equation ([5]) are determined by 



= 



(9) 



on dU. The evolution is determined by initial conditions 



n(x, 0) = no(2;), 0) = t>o(x), w{x,Q) = wq{x) 



(10) 



vq{x) and wq{x) are assumed to be nonnegative. Positivity and the integral of v and 
w are conserved throughout the dynamic. 

The equations of electro-hydrodynamics are the force balance equation of an incom- 
pressible Newtonian fluid coupled with a parabolic system of conservation equations 
and an elliptic equation. ([T]) and ([2]) are the Navier-Stokes equations with Lorentz force 
A(j)V(j). dni) and ^ are the conservation of mass equations of charge densities, u is the 
velocity field of the fluid, p is the pressure and v and w are densities proportional to the 
probability density of a system of negatively and positively respectively, self repelling 
charged particles in the fluid and (j) is the electrostatic potential due to Coulombic in- 
teractions of the charged particles along with the charge exterior to the domain fixing 
the boundary condition of (f). The electrostatic stress exerted by the charges on the 
fluid has the form of a rank one tensor and a pressure; 



Te stems from the balance of kinetic energy with electrostatic energy via the least action 
principle, |17j . For simplicity, we have assumed that the density, viscosity, mobility 
and dielectric constant are unity. 

The equations of electro-hydrodynamics are one of many fluid-particle systems 
which have attracted much attention for there challenge with regard to mathematical 
existence theory, derivation and simulation. Electro-hydrodnamic fluids are a particu- 
larly attractive complex fluid due to there emerging application in microfluidic devices, 
electric biogels, switchable soft matter and proton exchange membranes, [3l \T6\ [T8]. 
Some related works centered around other nonlinear Fokker-Plank and Navier Stokes 
systems and the equations of viscoelastic fluids can be found in [U [3 [HI [9l |T5] . In 
the case of electro- hydrodynamics, the closure of the nonlinear Fokker-Plank equations 
gives the Debye-Hiickel system, a basic model for the diffusion of ions in an electrolyte 
filling all of first formulated by W. Nernst and M. Plank at the end of the nineteenth 
century, p^. 

Main Result 

For the remainder of the paper, U is assumed to a connected, bounded open subset of 
with class C^'^ boundary dU. The space-time cylinder is Q = U x M"*". 
The first result of this paper concerns the existence of global in time solutions. 

Theorem 1. Given uq G H{U),vo G LP'{U),wq G L'^{U) there is a unique, global in 
time, weak solution u,v,w of ([I|)-([9]). 
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The definition of a global in time, weak solution will be given in the proof of 
theorem [1] at the end of section [2j The proof of theorem [1] is based on a modified 
Galerkin procedure found in [Rj. Formally setting u = 0, reduces (fUllOp to the Debye- 
Hiickel system. [5J have shown that the Debye-Hiickel system has a unique, global in 
time, weak solution under the assumptions on U above. It is not known whether global 
in time, weak solutions of the Debye-Hiickel system exist on general smooth, bounded 
domains in dimensions greater than 2. 

The second result of this paper concerns the long term behavior of the solutions 
guaranteed by theorem [TJ 

Theorem 2. There exist a positive constant A depending only on U and constant C| 
depending only on uq,vq and wq such that 

\Ht)\\H + \Ht) - VWh^u) + Mi) - W\\Uu) + U - nmiu) < Cf^~^'- (11) 

for all t £ where V, W, ^ is the unique steady-state solutions. 

The definition of the steady state solution and stationary solution are given in sec- 
tion [3l [5] have shown that solutions of the Debye-Hiickel system converge to the steady 
state solution in the norm in two and three dimensions provided the norm of solu- 
tions is bounded in time. Later, ^ proved that solutions of the Debye-Hiickel system 
on uniformly convex domains in arbitrary dimensions converge exponentially to the 
steady state solution, provided the solutions are defined globally in time. With slight 
modifications, the proof of theorem [2] presented in this paper implies the exponential 
convergence of solutions of the Debye-Hiickel system to the steady state solution in 
two and three dimensions without the assumption of uniform convexity on the domain. 
See corollary [1] at the end of section [3l 

Remark 1 (Some generalizations). In the forthcoming theory, the Laplacian in ([5]) may 
be replaced by any operator associated with a convex quadratic form Q{-,-) which 
is coercive over Hq{U) and satisfies estimates ()12p and (??) when is replaced 
by . In general, the system may include several density functions vi, . . . ,Vs with 
valences fii, . . . ,^s- Associated with such a parabolic- elliptic system is the entropy 



S{vi,. . . ,Vs) = / log + log Us 
and the interaction energy 

■^q{vi,. . . ,Vs) = min {Q{4>,^)- 4>{fiivi -\ \-fisVs)dx]. 



The kinetic energy is 



(S{u) = / \\u'^dx. 
u 2 



The entropy functional is 

J^q{u,vi,. ..,Vs) = S'{vi, ...,Vs)- ^q{vi, ...,Vs) + ^{u). 

If u is the motion of an incompressible, Newtonian fluid coupled with charge densities 
vi . . .Vs with interaction energy then J^Q(u{t),vi{t), . . . ,Vs(t)) is nonincreasing in 
t. In this paper s = 2, Q{v, v) = | Vup/2, //i = —1, vi = v, fi2 = 1, f2 = w. 
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Notation 



The LP{U) norm will be denoted || • ||lp{(7) while the Sobolev norm of W^''^{U) will be 
denoted || • ||i^i((7). The norm of a vector in Euclidean space will be denoted by | • |. 
We say a measurable function / is nonnegative if /(x) > for a.e. x G U. \\ ■ \\L\ogL{u) 
will denote the integral of / log / for a measurable, nonneagative function / provided 
that integral exists. See Chapter 3 of [19] for the defintion and properties of Banach 
space Xq valued functions / in U'{E; X) with derivative /' in L'^{E; X') for E an open 
subset of M and 1 < p,q < oo. If /' G L'^{E; X') then we take / to be the continuous 
representative of its equivalence class. 

In certain instances, it will be convenient to write the sum of the norm of two 
functions as follows; 

\\'^^v\Yx = \\u\\x + hWx 

with r = p when X = U'{U),H^{U) or U'{E;X) and r = 1 when X = L^{U) or 
LlogL{U). 

Let A-i : L'^{U) H^{U) n H^{U) be defined by A~^v = (f) provided = v and 
<t>eHl{U). 

A note on constants: we will use (const) to denote an inessential constant which 
may change from line to line. The letter C with various sub and superscripts will 
denote a constant refered to in various parts of the paper while constants ci, C2, C3, ca 
or defined below will be written to indicate which inequality was used in that line. 

Constants ci,C2 and C3 appear in the following versions of the Nash, Poincare and 
Sobolev inequalities resp. in two dimensions (see [13j): 

Iblli2([7) < ci\\v\\hi{u)\\v\\l^(^u)^ G C~(C/), 
\\v\\l^u)<C2\\Vv\\l,^^y yvGC^iU), 

\\v\\l3{U) < C3|bllJ/f(,7)lblli2^([/), G C~(f/). 

Constants ca and stemming from regularity of solutions of the Poisson equation on 
domains with C^'^ boundary (see [11]) will also be useful 



\\A-^v\\h^u) < ca\\v\\l^^u), Vt; G C^iU), 

(12) 

\\VA-^v\\Le<cML^\u)\\'^^~'''\\LHur V^;GC7~(^). 
Finally, for e > 0, [5j have shown that there exists > depending only on U 

M\h <4v\\ h^u)M\l log L{U)+Ce\\v\\ Li (U), G C°°(C/). (13) 

provided U is an open subset of with C^'^ boundary. 

The space of smooth, compactly supported, divergence free vector fields is 

r{U) = {ve {C^{U)f : V • = 0} . 

H(C/) is the completion of 'f{U) is the completion of Y{U) in the L^-topology and 
y{U) is the completion of Y{U) in the topology. {Ui}'^i is the L^-orthonormal 
basis of H(?7) of the first component of eigenfunctions of the Stokes operator with 
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eigenvalues {Aj}^]^. The dual space of V(C/) is V'(C/). See Chapter 1 of [T9] concerning 
these definition. 

The constants /i^, Mq, i?0; -^o '^o, which depend only on uq,vo and wq will 
later be important; 

fiv = vo dx, 
Ju 

/J'w = Wo dx, 
Ju 

Mo = ||uo||h(c/), 

Ro = max|||t!o||i2(f/), ||w^o|li2(c/)} , 
R'q = max{||?;o||LiogL(i/)) ll'"^o||LiogL(;7)} i 
= \\VA'^{vo - wo)\\l2i^u) 

are assumed to be finite. 



2 Weak Solutions 

In this section is a positive integer and to is a positive real number. 

Bi C L'^{[0,to];L'^{U)) is the ball of radius R and B2 C C{[0,to];W) the ball of 
radius M in their respective topologies. 

Define i : C{[0,to];W) ^ C{[0,to];r{U)) by 

V 

\{a)=^aiVi (14) 

and define ] : B\ ^ B\ ^ B\ ^ B\ ^ B\hy 

](v,vS) = {v,w,\l A'^iv -w)). (15) 

Note that i is an isometry from C([0, to]; IR'') into C([0,to]; H([/)). 
Define two operators 

: Si X Si X Si ^ ^2 
^ : S2 ^ Si X Si 

as follows; 3^(v,w,e) = a provided a is a solution to the dimensional system of 
ordinary differential equations 

a[ + AjCj + ajttf^ / Uj ■ VUk ■ Uidx = / {v — w)e ■ VUi dx (16) 

ai(0) = / uoUi dx, for a\\ i = 1, . . . ,u and t G [0, to] (17) 
Ju 
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and IV (a) = {v^w) provided {v^vS) is a weak solution of 

dv 

+ (i(a)) . Vv = V • (Vu - uV(/)) 

+ (i(a)) • Vi(; = V • (Vw + wV<j)) 

A0 = f — for t G [0, toll 
v{Q,-)=vo{-), w{f),-)=wo{-). 

Finally define 

Let be the constant (which depends only on U) specified in lemma [2] and [0,t') 
be the interval of existence of the equation /' = C^f"^, /(O) = max{i?Oi !}• Let Cq be 
the constant (which depends only on Rq and C*) specified in lemma [21 Let R = C^^^. 
FixO<t* <t' and let 

= (Mo + i?^)e**. 

For these choices of M, R and t* , which depend only on the initial data and the domain, 
it will be shown that 2f has a fixed point ai, when to = t* ■ Then, it will be shown 
that there is a constant Mi depending only on the initial data and the domain, but 
independent of M,R and to such that any functions a^,Vi, and Wu corresponding to a 
fixed point of ^ are bounded in the l'^ and norms resp. 

N.B. in the proof of proposition [3] we will make a slight abuse of notation by 
assuming uo,vq,wo used in the definitions above are not necessarily the same as those 
functions given in the introduction, resulting in possibly different M,R and to- 

Lemma 1. ^ is well defined, continuous and ^(Bi) is precompact in C([0, to]; M'^). 

Proof. (Well Definedness) Let (t;,t(;,e) E Bi. The system ([T6|), pT|) is a finite dimen- 
sional system of ordinary differential equations with continuous dependence on t and 
a. By Peano's theorem, there exists e > so that a solves (fTH]) . (fTT]) for t G [0, e]. We 
no show that a extends to [0, to]. Multiplying (I16p by Oj and summing over i = 1, . . . ,v 
we find that 

yi + -^i^i = yi^i [ iv-w)e- Ui dx 

< ( / \{v — w)e\'^ dx^ I ajUil'^ dx 

<\j^{v-w)Ux + \j^\e\Ux + \± 



1/2 



a?. 



Letting oj{t) = Y^^i=i ^li^) see that u> satisfies the differential inequality uj' < f{t)+uj 
with w{0) < Mq and Jq" f{s) ds < R^. We infer that uj{t) is majorized by the function 
(Mo + i?^)e* on the interval [0,to] so that 

V 

^ a\{t) < (Mq + i?^)e*o = M^, for ah t G [0, to]. 

1=1 
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Hence a G i?2 for all t = [0, to]- 

(Continuity) Since J^is well defined, a,a £ B2 provided {v,w,e),{v,id,e) £ Bi. 
Thus |aj(t)|, |aj(t)| < \/M for all t G [0, to] a-nd alH = 1, . . . , z^. Let 

= max max |?7j • V?7,- ■ Ud ■ 

i,j,k=l,...,u xeC\os{U) 

Subtracting the equations solved by a and a resp. from each other, multiplying by 
— and summing over i = 1, ... ,1/ we find that 

i=l 

» 1/ 2 

= / ^ ^ (ojOj - diaj){ak - aujUi ■ VUj ■ Uu dx 
+ / {{v — w)e — {v — w)e) {ai — ai)Uidx 

V 

i=l 

^ \J ~ ^-'^ ~ ~ '^^J [J ~ di)Ui\ dx j 

< (const) y^{(ii — di)^ H — / |(f — w)e — {v — ?D)ep dx. 

^=l ^ 

Letting r/(t) = ^J'=i(aj(t) — aj(t))^ we see that r/ satisfies the differential inequal- 
ity 77' < (conts)?7 + g{t) with r/(0) = 0. We have then for all t £ [0,to] \vi't)\ < 
g(const)to jto i^^g^i where Jq" \9{s) \ ds converges to zero as v, w, e converges to v, w, e 
in L4([0,to];^^). 

( Compactness) Fix i £ {1, . . . ,1^} and integrate (fTHl) over [s, t] C [0, to] to find that 
\aiis) - ai{t)\ 

= / / 2, (^j{''^)0'k{f)Uj ■ ■ Ui + {v — w)e ■ Uidx — Xiai{r) dr 
< {u^M^K^ + max XiM){t - s) + R^^/t - s < K'^/t - s 

i=l,...,u 

where K' = {v'^M'^Ky + v[iixx.i=\^ u \j^lvl^ -\- Flj^ provided |t — s\ ^ 1. For 6 > 0, choosing 
5 = min{l, e^/ {vK')} will show that \a{s) — a(t)| < e for all \t — s\ < 6 and all v, w, e in 
Bi . Consequently the image of Bi is uniformly equicontinuous in B2 . The compactness 
asserted by the lemma follows from the Arzela - Ascoli theorem. □ 

Lemma 2. '3^ is well defined and continuous. 
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Proof. Let a e B2 and write u = i(a). Let N e Z'^ to be chosen below and h = to/N. 
Define a sequence of functions in H^{U) as follows; for i = Z+ define = u\t=ih let 
and solve 

v'= (l + hA + hu' -V -hV ■ (v'V6'-^))~\'-^ 



= {l + hA + hu' -V -hV ■ {v'V(l)'-^)y 

= {I + hA + hu' -V + hV ■ {w'V(t>'~'^)y 
A(f)i = Vi-Wi 



with and for z = let = vojiv*^ = wq. The sequence is well defined by the Lax- 
Milgram theorem for example, and in fact satisfies 

||vi,u;i||2 + h\\Dvi,Dwi\\l 

= [ v'-^v' + w'-^w' dx + h f Vv' ■ V4>'-^Vv' - Vw' ■ V4>'-^Vw' dx 
Ju Ju 

=1 + II. 
We have by Holder's inequality 

i<lh\w'\\l + hv'-\w'-^\\l 
- 2" " 2" " 

and 

<!^\\Vv\Vw'\\l+^^{\\v'hM\v'\\2 + \\w%,2W 

< — ||Vi;\Vu;^||2 + /i(const)||t;\wii(||i;i_i,u;i_i||^ + 

If = ||t'*,ui*||2 and fii = ||Vu*, Viw^llg then 

LOi + h/ii < LOi-i + /i (const )u;j(u;j_i + (18) 

Let C* equal 2 times the constant in the previous line and / solve /' = C*/^ with 
/(O) = max{i?o, 1}. Then, for all t G [0,to], f{t) > 1, f'{t) is nondecreasing and 



1 f^'^'^ 1 

(const)/(t)(/(i) + fit)) < a fit) = fit) <-J^ f'is) ds = -if it + h)- fit)). 



rt+h 
It 

Hence 



fiih) + /t(const)/((z + l)h)ifiih) + fiih)) < fiii + l)h), 
for alH = 0, . . . , AT — 1 for /t sufficiently small and so the inequality 

Ui < coi-i + /i(eonst)a;j(a;j_i + ujf^i) 

< fiih) + /i(const)/((i + l)th)ifiih) + fiih)) < fiit + l)h) 
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for all i = 0, . . . , follows by induction on i with the case i = being ujq < Rq < /(O) 
for sufficiently small h. Summing equations (jlSp for i = 1, . . . N then implies 



N 

UN 

i=l 



JL rto 
+ y^hfii<Ro + (const) / f{t){f{t) + f{tf) dt 



Define functions Vh,Wh G L'^{[0,tQ]; H^{U)) as follows; 

Vh{t) = v\ Wh{t) = w\ for ih<t <{i + l)h and i = 0, . . . , - 1. 

Then 

f \\vh,WhfHidt= V h{uj, + n,) < Co(to) := i?o + (const) f{t){f{t) + f{tf) dt. 
Jo -^0 

The right-hand side of this inequality is bounded independently of h. Consequently 
some subsequence {vh,Wh}hio converge to {v,w) in the weak topology of L^([0,to]; 
H'^iU)) and the weak-* topology L°°{[0,tQ]; L^{U)). {u,v) satisfy 

\\v,w\\LH[0,to],LHU)) < lb'^llL2aOA.];i?HC/))ll^'^llL-([0,to];L2(C/)) " '^O^^" ^^^^ 

Let ip e C°°([0, to] X U), with Supp(V') C [0, to) x U. Write V*(-) = V'(«^, •) and compute; 

N-l ^(i+l)h /. N-1 



/ / ^'vhdxdt^Y^ / / i^'v'dxdt^Y^ / + '^)h) - il){ih)y dx 

Jo Ju Jth Ju ,;^o 

= - / Voworfa^-V" / ip' {v' - Vi-i) dx 
Ju ■/[/ 



N /■ii+l)h 



JU Jih Ju 



-0* dx dt. 



The shift operator is continuous on L'^{[0,tQ]; H^{U)). Thus choosing possibly another 
subsequence of /i | such that (ph converges to (p in L^{[0,to]] J^^{U)), the equation 
converges to 

ip'v dxdt = / V(0) 2^)^o dx + {Vv — vVcp) ■ Vip — u ■ Vvtp dx. 
QtQ Ju JQto 

The analogous equation holds for w and Acp = v — w (or almost every t S [0, to]- 

The well definedness (uniqueness of v, w) and the continuity of "3^ now follows from 
standard estimates on the solutions v,w (see for example [5j.) □ 

Proposition 1. ^ has a unique fixed point a^. 
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Proof, to, M, Rave defined so that ^ {j{Bi x Bi)) C B2 and W{B2) C Bi x B^. Thus 
2^{B2) = J2r o j o '3^[B2) C i?2- Furthermore, is a composition of two continuous 
maps followed by a compact map and hence is compact. A fixed point ai, of ^ in B2 
is guaranteed by the Schauder fixed point theorem. The uniqueness of the fixed point 
is proved by standard estimates. □ 

Proposition 2. Let < io < 00 and suppose that {vy,Wi,) = '3''{a^), and = 
J^{j{v^,Wi,)) exist for some choice of M and R. Then there exists constant Mi de- 
pending only on Mq, Rq, R^, Sq, fj.^, fj,^ and U, hut independent of M,R and to such 
that 

max |a,,(t)p + ^^^^/|| r2^m < Mi. 
Proof. Let <pi, = {vi, — w,^) and Ui, = i{a,^). Then one readily checks that 

'-^ (20) 
" I ft f ^ ^ 

i=i ^ Jo Ju 

for all t € [0,to]- 

Let f{t) = \\{v^)-, (u^iy)- ||^2((y)- / is absolutely continuous with /(O) = and one 
readily checks that f'{t) < for t sufficiently small. Consequently f{t) = for all 
t G [0, to] and Vi,{t) and Wi,{t) are nonnegative for all t G [0, to]- 

For /i > let 



jeh{t) = [ Mt) + h) \og{v^{t) + h) + {w^{t) + h) logK(t) + h) dx 

Ju 

\\VMt)\'^ + Mt){Mi) - Mt))) dx. 



2' 

J^fh is similarly absolutely continuous and 

d 



Xhit)= / (log(?;^ + h) + l- (f)^)vl + (log(u;^ + h) + 1 + (f)^)wl dx 
Ju 

+ / V(?!)i, • V(f>l - (f>l{vi, - w„) dx 
Ju 

Note that the last three terms vanish for almost evert t G [0, to] because (f)' = A~^{vl — 
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w^) for almost every t G [0,to]. Continuing, 
d f 

-^'^h{t) = J (logivi, + h) + l- (j)i,)vl + (log(u;;, + h) + l + (t)u)w'u dx 

= - V(log(i;i. + h) + I - (j)^) ■ {Vv^ - VuV(t)y + VyUjj) dx 
Ju 

- / V(log(ii;,y + h) + 1 + (j),^) ■ {Vwu + WuV(l)u + WuUu) dx 
Ju 



u Vy + h 

— 3—- {Vwu + {wi, + h)V(f>u) {Vwu + WyV(j)u) dx 
u Wy + h 

- / {vu- Wy)V(j)y ■ Uy dx. 
JU 

Integrating this expression over [0, t\ for t < to and adding it to (j20p (the last terms in 
both expressions cancel) gives 

" 1 " 1 /■* /■ 1 

V-a2(f)+jr^(t)< V-a2(0) + jrft(0)- / / —\Vvy-VyV(l)y\^dxds 

7^ ^ ^ Jo Ju Vu + h 



i=l 



* r 1 2 

I Vtfjy + WjyV^jyl dxds + r](t) 



JU Wu + h 
where ^ 

'nit) =11 T^<t^u ■ (Vwj, - Vuy(t)v) dx dts 

Jo Ju Viy + i^ 

-V4>v iywu + WuVcfiv) dxds 

lo Ju + n 

is majorized by the two terms appearing directly to the left of it and 




JO Ju 



By bounded convergence we have 
" I " 1 

1=1 i=l 

<E^«?(0)+limJr,(0) + ^ f [ \VcP\^dxds (21) 
~l ^ MO ^ Jo Ju 
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for all t G [0,to]- 

We infer from (1211) that 



|a,(t)|2 + \K{t),w,{t)\\LiosLiU) + \\^Mt)\\hiu) < C := Mo + 2ii^ + 25o 
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for all t G [0,to]- We use this result to make an energy estimate of v^, and Wiy. We 
estimate 

v,yV(l)^ • Vvjydx < \\Vv,y\\L2(^u)\\vu\\L3{u)\\'^4>u\\L(i{u) 

u 

< c'^\\Vvu\\l2(u)\\vu,w„\\l3(^u)\\v^ -w„\\]^'i\\V(p\\]^2(^u) 

1 /2 

< c'^\/C'||Vt;i.||i2([7) (^e\\Vu,Wu\\H^^^\\v^,Wi,\\LiogL{U) + Ce\\Vu,Wy\\Li^u)^ 

< c'j^VC'\\Vvu\\l2(u) (c'€\\v^,w^\\]jj_f^,js^ + C^{fj.y + fi^ij 

I 

< 3c^C'^/e\\v,„W^\\j^^u^ + —VC'CeiHv + ^J'w) 

An analogous inequality holds for w^. Choose e sufficiently small so that 3c^C"-y/e < 1/4 
so that 

2 

2 

< 2ci\\Vv^,Vw^\\]^2(^u)\\v,y,W,y\\Ll(^u) + \\VVi,,VWi,\\l2(^u^ + -VC'CeilJ.v + fJ-w) 



LHU) 

3 2 



e 



\2 



Let M[ = 2{e)-^VC'C,{fi^+fi^)+2cl{fiy+tiu,)^ + 2Ro. Thenuit) = \\v^{t),w^{t)\\l2^^^ 
and 7] = \\\/viy(t),Vw,^(t)\\'j^2(^if^/2 satisfy 

—uj + oj + 'n<Mi, uj(0) < Ml 
at 

so that w{t) < Ml for ah t G [0, to]- Setting Mi = M[ + C the theorem is proved. □ 
Proposition 3. There exist functions 



such that for any < tQ < oo, {vy,Wy) = ^^(a,^), and a^ = ^ {j{vy,Wy)) and 

max (|a^(t)|2 + \\vu{t),w^{t)\\l2,u^ < Mi. 

Proof. Let T be the maximal positive real number for which the theorem holds. Clearly 
to > t* > 0. Then \au{to)\'^ ,\\vu{to),Wu{to)\\'j^2(^if^ < oo. Consider the functions a*, 
= i^ial) guaranteed by Theorem [1] when we introduce Uq = i(a*(to)), Vq = 
Vu{to), and Wq = w^{to) (also defining Mq, Rq, Rq and which are finite in terms of 
Uq,Vq and Wq.) Let ai,{t) = ay{t) for t < T and du{t) = al,{t — to) for t > T. Define Vi, 
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and Wy in the same way. Then there exists t* > with 0<to = ^ + ^* <oo and some 
M and R for which {vy,Wy) = '3^ {ay), = St^ {j{vy,Wi,)). By Theorem [2l 

max \K{t)\^ + \\vu{t),Wy{t)\\\2(u\ < ^h- 
te[o,to] ^ ' 

Then > T contradicts the maximaUty of T and concludes the proof. □ 
With the help of proposition [3] we are now able to supply the 



Proof of Theoremm {u,v,w) is said to be a global in time, weak solution of (fTl- llOp 
provided 

u G L2(M+; V(C/)) n C(IR+; V(?7)) with u' £ L2(R+; V'(C/)) 
v,we Ll^{R+;H\U))nC{R+;L'^{U)) with v',w' G lI^{R+; H-\U)) 

and if for all t G M+ and ah summable test functions / G C°°(M+; ^(C/)), g,h,ip G 
C7~(Clos(Q)) with ^\dU = 0, 

/ {u'J)dt+ u-Vu- fdxdt = - Vu-Vf -{v-w)V4>- fdxdt 
Jm.+ Jq Jq 

/ {v',g)dt+ / u-Vvgdxdt = - {Vv - vV4>) -Vgdxdt 
Jr+ Jq Jq 

/ {w',h)dt+ / u-\/whdxdt = - {Vw + wVc/)) ■ Vhdx dt 
Jr+ Jq Jq 

V(p ■ Vtp dxdt = — / {v — w)il) dx dt 



Q 

Let ay,Vy and Wy be as in proposition [3l Let T > 0. Let Uy = i(a). Then there is a 
constant Ct < oo such that 



rT 

max |n^(OlH(c/) + i^^^ WIIl2{(7) + \K{i)\\\/(u) + \\K{t)\\\/'{u) dt 



T 

2 



(22) 



We have then for some subsequence reindexed by ly 

u^^umL'^{[0,T];y{U)), n in L~([0, T]; H(C/)) (23) 

<-n'inL2([0,T];V'(C/)) (24) 

7;^,w;^^i;,u;inL2([0,T];Fi(f/)), t;^, u;^ i;, u; in L°°([0, T]; ^^(c/)) (25) 

v'y,w'y^v', w' in L2([0, T];H-\U)), (26) 

for some functions u, v and tf. We will see later that u, v and w do not depend on the 
T. By the Aubin-Lion compactness lemma (Theorem 2.3 in [l_9j), we have 

Uy^umL'^{[0,T];H{U)), Vy,Wy ^ v,w in L'^{[0,T]; L'^{U)). (27) 
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It follows from elliptic regularity that there is also (j) with (after choosing possibly 
another subsequence) 

(l),^(l,\nL\[Q,T]-H\U)). (28) 

Let g G C°°{C\os{Q)) with Supp((7) C [0,T) x U. By the triangular ineqaulity, the 
equality '3^{ay) = {vy,Wi,) implies 

/ {v',g)dt+ / u-Vvg+ / {Vv - vVcj)) - Vgdxdt = 
Jm.+ Jq Jq 

/ {vl,g) dt+ / ■ Vv^g + / {Vv„ - Vi,V(p,y) • Vg dx dt 
Jr+ Jq Jq 

/ {v',g)dt+ / u-\7vg+ / {Vv - vV<f>) ■ Vg dx dt 
Jr+ Jq Jq 



lim 







provided we can demonstrate the convergence of the individual terms. The conver- 
gence of the linear terms follows from the definition of weak convergence. We check 
convergence of the quadratic terms; 



lim 

i^— >oo 



/ {u,y ■ Vvu — u ■ Vv) g dx dt 

Jq 



- iH^"^'' ~ ^llL2{[0,T];H((/))l|VVi.||L2{[0,T];Hi(C/))SUp|fi'| 

Q 



+ lim 

l/—^00 



Q 



u ■ V{vp — v)g dx dt 



0. 



where the first term converges by (|22p . (j25p and (j27p . Similarly 



lim 



/ {vi/\/(j)u — vV4>) ■ Vg dx dt 

Jq 

- i™o "^'^ ~ '"\\L2{[0,T]-L^U))\\'^4>u\\L2{[0,T];m{U)) SUp|V5r| 

Q 

+ Jhn^ ll'^llL2([0,r];L2(l/))l|V((/'i. - (p)\\L2([0,T];m{U)) SUp |Vfif| = 0. 

The analogous limits hold for u>^ and h G C°°(Clos(Q)). Let ip e C°°(Clos(Q)) with 
tpldU = 0. By dSZl) and (f28l). 



lim / V(j)u • VV' da; c^i 



{v — w)il) dx dt. 

Jq 



Let F be of the form F{t,x) = J2i=j bj{t)U{x) for {bj{-)}j^^ continuous. By the 
triangular ineqaulity, the equality ^ {j{vy,Wi,)) = ay implies 

/ {u',F)dt+ / U-VU-F+ / Vu-VF - {v - w)V(j)- Fdxdt = 
Jr+ Jq Jq 

/ «, F)dt+ / Uj, • Vn^ ■ F+ Vuy ■ VF - (v^ - Wu)V<i)u -Fdxdt 
Jr+ Jq Jq 

/ {u',F)dt+ / U-VU-F+ / Vu-VF - {v -w)V(l)- Fdxdt =0 
Jr+ Jq Jq 



lim 

I/— >oo 
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provided we can demonstrate the convergence of the individual terms. The convergence 
of the hnear terms follows from the definition of weak convergence. Convergence of the 
first quadratic term can be found in We check the last term 



lim 



{{vu — Wy)V(t)y — {v — w)V(t)) ■ F dxdt 



- }^ "^'^ -V,Wu - w||L2([0,r];L2(;7))l|V</'i.||L2([0,T];L2(C/)) SUp|F| 
+ i™oll^ ~ W'I|l2([o,T];L2(;7))||V((?;»,, - (/')||l2([o,t];L2([/)) sup |F| = 0. 

Q 



by ([271) and (|28D . 
Finally, 

/ {u,f)dt+ [ u-Vu- f + Vu-Vf - iv-w)V(j)- fdxdt = Q 

for / E C°°([0, T]; Y) by approximating / by functions of the form Yli=j bjit)Ui in the 
C'^ topology for = 1, 2, . . . 

The uniqueness of n, v, w, which implies the extension of u,v,w to t ^ M+, and the 
fact that 

limii(t) = no, limv{t) = vq, limw{t) = wq, 



now follow from standard estimates. 



□ 



3 Long Term Behavior Behavior 

Define a functional ^ on Hq{U) and two absolutely continuous functions and ^ 
over R_|_ in terms of functions u, v, w, (j); 

= J ^\V4>\'^ dx + ^Xv^og ^ J ex.p(j)dx^ + fiy^ log ^ J {expcpy^ dx^ (29) 

^(t)= / v{t)logv{t)+w{t)logw{t) + hvcj){t)\^ + hu{t)\^dx (30) 
Ju 2 2 

^(0 = y^^i^i +2 V ^2 V ^\^m)-n^^^ (31) 

^ is the entropy function of electro-hydrodynamics while ££ is stems form a function 
introduced in [1] to study convergence of the Debye system to the steady state solu- 
tion assuming Dirichelet boundary conditions. is a positive constant which will be 
specified in lemma [6l ^(0 is strictly convex and bounded from below. There exists a 
unique function $ G C°^{U) n Co(C/) such that ^($) < ^(0) for ah $ / G UliJJ), 
e.g. |12j : Define functions 

^Ti \ exp$(x) exp(-<I>(x)) 

Jj^exp$(x)dx Jj^exp(-$(x)) 
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We call V, W, $ the steady state solution. We will frequently us the fact that there are 
constants a, b, a' , b' for which 

< a < V{x) <b<oo, 0<a' < W{x) <b' <co, Vx G U. 

We call V, W, <I>, U the stationary solution when [7 = 0. An important fact about ^ is 
that 

A^V^ = V{V + W) (32) 

is the gradient of a pressure so that the stationary equations are consistent. 

We first recall a well known fact about entropy functions which holds additionally 
with the kinetic energy term |tip/2 found in J^. The manipulations in differentiating 

may be justified by approximating v, w by strictly positive functions or by the 
argument used in the proof of proposition [2j 

Lemma 3. 



J^{t)-Jf {{)) = - I \2V^/v{s) - ^/v(s)Vcj){s)\'^dxds 
Jo Ju 



t 

Ju 



2Vy/w{s) + x/u;(s)V0(s)|2 + \Vu{s)\'^ dx ds < 0. 



Proof. The proof is contained in the proof of proposition [2j □ 

Lemma 4 (Weighted Poincare Inequality) . Let Q be an open subset of M"' and p G 
H\^}) satisfy 

< p{x) <b < oo, \fx ei} 
with p^^ integrable. Then there exists aCp = C{p,Q) such that 

f^dx<C^ [ |V(/p)p dx 
Jn 

whenever 

fdx = 0. 



Proof. Suppose that no such constant exists. Then there is a sequence of functions 
in H\n) with 

fidx = 0, \\fi\\ly^.>i [ \V{f,p)\^dx. 



lL2(f7) 

Let hi = fi/\\fi\\L2{n) so that 



lL2(n) 

Write gi = hi p. Then 



1 = > i I \Vih,p)\'^dx. (33) 



gf dx = hfp'^ dx < b^ 
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shows that gi is bounded in H^{Q) and thus converges weakly to an element g £ H^{Q). 
Fatou's lemma 

/ iV^l^dx < liminf / \Vgi\^dx< lim - = 
shows that g{x) = G a.e. for some constant G. Then 

G [ p^^ dx ~ 1 p^^g dx = lim / P^^gi dx — lim / hidx — lim -rr-rj, / fidx = 0. 

We infer G = since is nonzero on a set of positive measure. The contradiction 
with (j33p completes the proof. □ 

Lemma 5. Define functions g,h > via the formula 

g{x, t)V{x, t) = v{x, t), h{x, t)W{x, t) = w{x, t) 
for {x,t) gU X [0, oo). Then 

gV dx = fly, / hWdx = fi^, Vt G [0, oo) (34) 
u Ju 

and g and h are generalized solutions of the equations 

diaV) do d"^ 

-^l^+u.V{gV)=V-{V{Vg-gV^)), / - 5 = 0, (35) 

Ot Of Of 

^ ' ' + u- V{hW) = V • (ly (Vfe + W^)) , ^ + /i^ = (36) 
ot ov ov 

where we define ^ = (j) — ^. 
Proof. By definition 



gVdx = / vdx = fiy, / hWdx = / w dx = fiw 
u Ju Ju Ju 



Elementary manipulations of the definitions will show (j35l I36p . For simplicity assume 
g,v > are smooth (the alternate case being treated by approximation); 

log f = log 5 + <I> + log f , = log 5 + <I> + (const) . 

Thus, 

-Vg - = -Vv - V$ - = -Vv - \7(t>. 

g V V 

This shows the second equality in (135]) . Multiplying this equation hy v = gV, and 
taking the divergence 

V-{V{Vg-gV^)) = ^^+u-Vv = ^^ + u-V{gV) 
as required. The analogous calculation proves the result for h. □ 
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Lemma 6. There exists positive constants Ci,C2 and C3 depending only on U, V, 
ll^llinfi ll^llsup, W, llW^llinf and ||VF||sup such that 

j^^it) < -Ci^it) + C2^{tf + C^J^{t)\ (37) 

Proof. Without loss of generality, we may assume w = W = the general case simply 
being a sum of the argument given below for v,V to that oiw,W. Define E = v — V and 
^ = (j) — ^. Using elementary manipulations and noting that E = V{g — {E /V) = 
Vg, dtE = Vdtg we find 



- ^V^ • - 2V^ • V^- + 2 ("^ + 1 ] IV^-I^ 1 dx 



+ [ -G| Vup + Ou ■ V0A(/) + ■ V{E + V) dx. 
Ju ^ 



lU 

By application of Young's inequality, the last three terms in the first integral may be 
bounded in terms of 1/2 the first term plus products of terms involving \E\ and |V^'| 
of order greater than 2. Since A$V$ is a gradient, (j32p . we may insert the term 
9m • V^A^ in the second integral. Then u ■ V^A^f = u ■ VipAcf) - u ■ V^Acf) - u ■ 
V0A<I> + u ■ V<I>A$ shows that 



u-V(pA(pdx = / n- V^A^ + u- V$A0 + n- V0A$-n- V$A$dx 
Ju 

= / u-V'i/A'^ -u-{A^V(j))+u-V(j)A^dx< / ai\u\^ + a2\E\^\V'^\^ dx 
Ju Ju 

where 4aia2 = 1. We may estimate the third term in the second integral by rewriting 
it as 



E f 
-u-V-E + u-V^Edx < / hlul'^ + bilEl^ + ibs + Cu.h) 
V Ju 



E 



+ h\\^^\\sup\u\'^ dx 



where 46162^3 = 1 and 46465 = 1 and C^j is the constant in lemma|3]with in place 
of p and E (which has integral zero) in place of /. Choose Qai, 63 and 65 sufficiently 
small so that 

C4 := min |^||y||min - 63 - Cujh, 1 - 2eaiC2| > 0. 

With 61, . . . , 65, oi, 02 specified we fix = 64II V<^>||g^p. Applying lemma H] and the 
Poincare inequalities once more implies 



|^(t) < -C'.imt) + ||Vm||2 + \\V{E/V)\\U 

+ C5 f lul^ + lEl^ + E^lV^"^ + \E\\V^'^dx 



where C4, C5 = C4, (75(61, . . . , 65, ai, 02,02, Cuj, C4, ||l^||sup, ll^llinf)- [T] have treated the 
Dirichelet boundary condition case to arrive at an inequality similar to this one, but 
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the Gahardo-Nirenberg-Sobolev inequahty and regularity of solutions to the Poisson 
equation are sufficient to arrive at the following; using the estimates in the proof of 
theorem 1.3 in [1] one may show that there are constants C(, C2 and C3 so that 

C5 [ \u\^ + \E\^ + E^lV^l"^ + lEWV^l"^ dx < 
Ju 

C[e[J^{t) + \\Vu\\l + \\V{E/V)\\l,] + ^^{tf + ^L{tf 

for all e > 0. Choosing eCi < C4 completes the proof. □ 
Lemma 7. There exists a subsequence ^1"^+ f^''' which 

lim \\v{tj) - V\\l2 + +\\w{tj) - W\\l2 + ||V(/.(t,-) - ^'11^2 + ||n(t,-)||i2 = 0. 

Proof. [5] showed in the proof of theorem 6 of that article that lemma [3] is sufficient 
to find a subsequence of M+ for which the first three limits hold. One must 

assume that supjg^ tt;||/^2 < co which is the case here. Since niay be cho- 

sen from a set of positive measure, by lemma [3l we may choose a second subsequence 
{*j"}j"=i {^j'}j?=i foi' which additionally limj//_^oo l|V^(*i")llL2 = 0. By the Kon- 
drakov's embedding theorem (||ii(i)||L2 is bounded by K{{}) < 00 for all t G R_|_) there 
is a third subsequence {tj}^^ for which u{tj) converges to function [/ E V. By Fatou's 
lemma, ||VC/||^2 = so that [/ is a constant function. But C/ € V then implies that 
U = 0, giving the fourth limit. □ 

Lemma [6] and lemma [7] combined imply the exponential convergence to the station- 
ary solution. 

Proof of theorem\^ Let be the sequence provided in lemma [71 It follows from 

lemma [6] that \\m.j^ooL{tj) = 0. Choose j sufficiently large so that 



C2L{tjf + C^L{tjf<-^L{t,). 



Then for t >tj, L{t) is nonincreasing and so 

j^L{t) < -CiL{t) + C2L{tf + C^L{tf < -^m. 

Setting A = Ci/2 and Q = max {l, HFUj^jj^, HVFllj^i^} sup^gfo,*,] -^(t) completes the 
proof. □ 

The proofs of lemmas [3]through[7] and the proof of theorem[2]hold similarly solutions 
of the Debye-Hiickel system by formally setting n = 0. We thus have 

Corollary 1. There exists a positive constant A depending only on U, V, ||l^||inf, 
IIV^II supi ||VF||inf cind ||VF||sup such that the solution to the Dehye-Huckel system 
V, w, (f> in dimensions 2 and 3 converges with rate e~'^* to steady state solution V, W, $ 
in the L^,L^ and norms respectively provided supjg ||f,ti;||/^2 < 00. 
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Remark 2. This corollary in some sense improves theorem 2 of JJI where the domain 
U is assumed to be uniformly convex. IfUis uniformly convex, then one may produce 
a bound 

i^reiW := [ v\ogv + w\ogw-V\ogV -WlogW + ^\V^''dx>-^j^K,,,{t) 

for some X' depending on U by applying remark 3. 7 of 12] concerning a logarithmic 
Sobolev inequality for bounded domains and noting the convexity of J. The conver- 
gence of V, w to V,W follows then from the Csisdr-Kullback inequality, 120^ . 

4 Conclusion 

The equations of a incompressible, Newtonian fluid coupled with charges in two dimen- 
sions have been studied. The key step toward the existence of global in time solutions 
is the existence of a decaying entropy function which guarantees the dissipation of ki- 
netic and electrostatic energy and entropy. Future avenues of study are the regularity 
of these solutions and generalizations to other incompressible Newtonian fluid systems 
coupled with polarized particles. The state of knowlege (rather lack of knowledge) 
concerning the global in time existence of weak solutions to the Debye-Hiickel system 
in three dimensions prevents the generalization of theorem [1] in this paper to three 
dimensions. 
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